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Abstract 

Using 20.7 pb~^ of e+e~ annihilation data taken at s/s = 3.671 GeV with the CLEO-c de- 
tector, precision measurements of the electromagnetic form factors of the charged pion, charged 
kaon, and proton have been made for timelike momentum transfer of = 13.48 GeV^ by the 
reaction e^e~ h'^h^ . The measurements are the first ever with identified pions and kaons 
of |Q2| > 4 GeV2, with the results F^(13.48 GeV^) = 0.075 ± 0.008(stat) ± 0.005(syst) and 
Fft:(13.48 GeV^) = 0.063±0.004(stat)±0.001(syst). The result for the proton, assuming = G^, 
is G^(13.48 GeV^) = 0.014 ± 0.002(stat) ± O.OOl(syst), which is in agreement with earlier results. 
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Electromagnetic form factors of hadrons are among the most important physical observ- 
ables. They provide direct insight into the distribution of charges, currents, color, and flavor 
in the hadron. Form factors for spacelike momentum transfers, > 0, are determined by 
elastic scattering of electrons from hadrons available as targets. Form factors for timelike 
momentum transfers, < 0, are measured by annihilation e~^e^ ^ h^h^. They lead to 
insight into the wave function of the hadron in terms of its partonic constituents. In this 
Letter we report on the first precision measurements for the timelike electromagnetic form 
factors of the pion, kaon, and proton, for = 13.48 GeV^, by means of the reactions 

e^e~ TT^TT^, K^K~, and pp (1) 

Measurements of the timelike form factors of pion and kaon, and Fx, with identified 
pions and kaons exist for \Q^\ < 4.5 GeV^ [7]. For larger either only upper limits for 
F{\Q^\) exist, or the few observed hadron pairs, not separately identified as pions or kaons, 
are divided between the two according to the expectations based on the vector dominance 
model (VDM) to obtain F^, and F^ 

Timelike form factors of the proton for \Q'^\ > 6 GeV^ were first measured by the Fermilab 
E760/E835 experiments via the reaction pp e^e~ According to perturbative QCD 
(pQCD) at large momentum transfers the timelike form factors of protons are expected 
to be nearly equal to the spacelike form factors, but the Fermilab measurements found 
the timelike form factors to be nearly twice as large. In this Letter we provide the first 
independent confirmation of the Fermilab observations. 

Theoretical predictions of form factors based on pQCD rely on the validity of factorization 
for sufficiently high momentum transfers, and lead to quark counting rules, which predict 
0] that F{\Q^\) (X IQ^I^"", where n is the number of quarks, so that F{\Q^\) oc as/\Q^\ for 
mesons, and oc a|/|<5''| for baryons. pQCD also predicts that the form factors 

for the helicity-zero mesons m = vr, K, p, ... are proportional to the squares of their decay 
constants so that F^{\Q^\) / FkHQ^]) = fl/fh as \Q'^\ ^ oo. 

By providing tests of the above predictions, we expect to shed light on the important 
question of the momentum transfers which are sufficiently large to validate the use of pQCD, 
a question which has been in debate for a long time . 

The timelike form factors of the charged helicity-zero mesons are related to the differential 
and total cross sections for their pair production by 

^(e+e- - m+m-) = ^^^1 |F^(s)| W^, (2) 

where m = vr or i^', a is the fine-structure constant, (3m is the meson velocity in the laboratory 
system, s = is the center-of-mass energy squared, is the electromagnetic form 

factor of the meson, and 6 is the laboratory angle between the meson and the positron beam. 

The e+e~ —>■ pp differential cross sections are related to the magnetic (G^) and electric 
(G^) form factors of the proton. With r = Arrip/ s 

^ = ^/5p[|G^M(^)P(l + cos^^)+r|G^^(.)Psin^^] (3) 

In the present measurements we do not have sufficient statistics to separate G^^{s) and 
G%j{s). We therefore analyze our data with the two assumptions, G^^ = G^j and G^^ = 0, 
as in Ref . 



3 




0.80 0.85 0.90 0.95 1.00 1.05 1.10 1.15 1.20 



X, = (E,+ + E,.)//i- 

FIG. 1: MC distributions as a function of = {Eh+ + Eh_)/^/s, where h = tt (top), h = K 
(middle), and h = p (bottom). The signal regions are defined as 0.98 < < 1-02 and are 
designated by the vertical arrows. 

The e~^e~ annihilation data samples used in the present measurements consist of 20.7 
pb-i at = 3.671 GeV and 2.89 pb^^ at the V(2^) mass, ^/s = 3.686 GeV. The data were 
collected at the Cornell Electron Storage Ring (CESR) with the detector in the CLEO- 
c configuration The detector is cylindrically symmetric and provides 93% coverage 
of solid angle for charged and neutral particle identification. The detector components 
important for this analysis are the wire vertex detector (ZD), the main drift chamber (DR), 
the ring-imaging Cherenkov detector (RICH), and the Csl crystal calorimeter (CC). They 
are operated within a 1.0 T magnetic field produced by a superconducting solenoid located 
directly outside of the CC. The ZD and DR detect charged particles, and the DR provides 
measurement of their momenta with a precision of ~0.6% at p = 1 GeV/c and ionization 
energy loss (dE/dx). The RICH detector provides particle identification, and covers 80% of 
the solid angle. The combination of dE/dx and information from the RICH detector allows 
for separating different charged particle species. The CC allows precision measurements of 
electromagnetic shower energy and position. 

A fully reconstructed event is required to have two charged particles and zero net charge. 
The charged particles in the Tr"*"??" final state analysis must have |cos6'| < 0.75 and have 
an associated shower in the CC. The charged particles in the K^K~ and pp analyses must 
have I cos^l < 0.80. Each charged particle is required to satisfy standard requirements for 
track quality and distance of closest approach to the interaction point. For the tt+tt" and 
pp analyses, the two charged particles must have a net momentum < 100 MeV/c, and for 
the K^K~ analysis, they must have a net momentum < 60 MeV/c. 
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FIG. 2: Data events as a function of for vr+vr (left), K^K (middle), and pp (right) final 
states. The dashed lines denote the signal regions defined as 0.98 < Xh < 1.02. 



The hadronic e^e^ h^h^ processes, where h = 7r,K,p, are studied by generating 
signal Monte Carlo (MC) samples, using GEANT P| to simulate the CLEO-c detector. 
MC samples of the leptonic background processes e+e~ {I = e, /i) are also studied. 

Figure H shows the MC distributions for two track final states which pass the tt+tt", K^K~ 
and pp criteria described above. The normalized center-of-mass energy variable is defined 
as the sum of the energy of the two tracks (assuming the particle species of interest for each 
track) divided by y/s. Figure ^ shows that the e^e~ — > K^K~ and e+e~ — >■ pp signal 
regions are sufficiently displaced from the dominant e"'"e^ l^l^ background, while the 
e^e^ 7r+7r~ signal overlaps with e"'"e^ l^l^ . 

In order to suppress background events, it is first required that the accepted events 
have the ratio of the CC determined energy Ecc and the track determined momentum p, 
Ecc/p, be less than 0.85. For pp events this cut is only applied on the positive track, for 
7r+7r~ and K^K~ it is applied on both tracks. In order to obtain a higher level of lepton 
rejection, signal to background optimization studies are made in terms of a likelihood variable 
defined by using RICH and dE/dx information, L{p, K) — L{1) = Lrich{p, K) — Lrich{1) + 



TABLE I: Summary of form factor results. The first errors are statistical only. The second errors 
in cross sections and form factors are systematic errors as discussed in the text. The form factor 
results for protons correspond to the assumption G^^ = G^^. Results for the assumption = 
are ~9% larger. 



pion kaon proton 

# of Observed Events 26.0 ± 5.1 72.0 ± 8.5 16.0 ± 4.8 
Lepton Contribution ~ 0.1 0.6 ± 0.2 < 0.1 
'0(25') Contribution < 0.1 0.6 ±0.1 1.9 ±0.2 

# of Net Signal Events 25.9 ±5.1 70.9 ± 8.5 14.1 ±4.8 
o-(e+e- ^ /i+/i-), pb 9.0 ±1.8 ±1.3 5.7 ±0.7 ±0.3 1.2 ±0.4 ±0.1 
Ffe(|Q2|) F^=0.075±0.008±0.005 Fr-=0.063±0.004±0.001 Ga/=0.014±0.002±0.001 
Q"F,,(|g2|)^ GeV^ |Q2|i?^=i.oi±0.11±0.07 |Q2|i?^=o.85±0.05±0.02 |Q^|GAf=2.53±0.36±0.11 
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^dE/dxiPy ^) ~ '^dE/dxi^)- ^he Optimum requirement is found to be L{p, K) — L(e) < 0, and 
L(p, K) — L{fi) < —2 for each track. 

Rejecting leptonic background in the tt+tt" sample requires additional measures. These 
are determined by studying radiative Bhabha events in the continuum {^/s = 3.671 GeV) 
data, by studying /i tracks from the e^e^ ^^^^ Monte Carlo sample, and by studying 
pions of appropriate momenta (~ 1.83 GeV/c) in the existing CLEO sample of inclusive 

— > K~TT^ data taken at ^/s = 10.58 GeV. The optimization criteria which emerged from 
these studies are that pions must deposit Ecc > 0.42 GeV, and must have L{n) — L{e) < 0, 
and L{K) - L(7r) > 0. 

Figure |21 shows the distribution of the events which meet all the above selection criteria 
as a function of X^. The signal region is defined as 0.98 < < 1.02 as bounded by the 
dashed lines. The observed counts, the estimated lepton contamination counts, the counts 
contributed by the tail of the ip{2S) resonance, and the net signal counts are listed in Table 

HI 

The net signal counts are related to the Born cross sections as ao{e~^e h^h ) = 
N/[eC{l + 6)], where e is the detection efficiency, C is the total integrated luminosity, 
and (1 + 6) is the radiative correction factor associated with h^h~ production from e^e~ 
annihilations. The proton and the kaon detection efficiencies, = 0.657 ± 0.003, and 
ex = 0.743 ± 0.003, are determined from the signal MC samples. The pion detection effi- 
ciency, e^r = 0.166 ± 0.002, is determined by signal MC sample and the —>■ K~tt^ data. 
The signal MC samples are generated with angular distributions according to Eqs. |21 and 
121 The initial state radiation corrections are determined using the method of Bonneau and 
Martin with the addition of /i and r pair loops to the vacuum polarization term. For 
the final states vr+vr", K~^K~ , and pp, the (1 + 6) correction factors are 0.832, 0.810, and 
0.860, respectively. 

The resulting cross sections are listed in Table HI Various sources of systematic uncer- 

TABLE II: Sources of systematic uncertainty for the e^e~ — > TT~^ir~ , e^e~ — > K^K^ , and e^e^ 
pp cross sections. 



Source 


TT+TT- (%) K 




pp (%) 


Trigger 


1.0 


1.0 


1.0 


Tracking 


2x1.0 


2x1.0 


2x1.0 


Xh Signal Region 


4.1 


0.5 


3.7 


Net Momentum 


4.8 


2.6 


6.9 


dE/dx+mCR PID 


2x2.7 


2x1.2 


2x1.6 


Ecc 


10.7 






en{Ecc) 


2x2.3 






MC statistics 


1.3 


0.4 


0.5 


tp{2S) Contamination 


0.08 


0.1 


1.0 


Leptonic Background 


0.05 


0.3 


0.0 


Radiative Correction 


0.2 


0.2 


0.2 


Luminosity, £ 


1.0 


1.0 


1.0 


Total (in quadrature) 


14.6 


4.4 


8.9 
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FIG. 3: Compilation of the existing experimental data for the pion (top), kaon (middle), and proton 
(bottom) form factors with timelike momentum transfer. Top and middle: The solid points are 
from identified vr^ and 0. The open points are from unidentified h^, divided into vr^ and 
according to VDM expectations. For pions (top) the open triangle denotes the value obtained at 
M'^^J/ip) in Ref. [l^]. Bottom: The solid points for proton form factors are from the measurements 
and compilation of Ref. 0|. The arbitrarily normalized solid curves show the variation of a5(|Q^|) 
(top and middle) and a^(|(5^|) (bottom), as determined for four flavors and A = 0.322 GeV. The 
dashed curve in the bottom plot shows a|(|(5^|) fit to the spacelike form factors of the proton. 



tainties in the cross sections have been studied. These are listed in Table |n] Their sum in 
quadrature is 14.6% for pions, 4.4% for kaons, and 8.9% for protons. Integrating Eqs. |21and 
El leads us to our final results for the form factors as listed in Table HI 0] . 

The 2.89 pb~^ of data taken at ^/s = 3.686 GeV, or the ip{2S) resonance, were used in 
the form factor analysis to evaluate the contribution that the resonance makes to the form 
factor data at ^/s = 3.671 GeV. For this purpose the ^'(25') data were analyzed in exactly 
the same way as the form factor data. Our yields in the ip{2S) data sample, which we 
use to estimate the background from the ip{2S) tail feeding into the continuum sample, is 
consistent with the expectation based on the branching fractions for p and K, but lower 
for TT. 

Our results for timelike form factors are displayed in Fig. |21 as \Q'^\Ft^, \Q'^\Fx, and 

together with the existing world data for the same. 
Our precision result for the pion, |Q2|i7;(i3.48 GeV^) = 1.01±0.11±0.07 GeV^, is the 
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first such directly measured result. It provides empirical validity for |Q^|F^(9.6 GeV^) = 
(0.94±0.06) GeV^ obtained by interpreting T{J/ijj —>■ 7i^TT~)/T{J/ip —>■ e^e~) as a mea- 
sure of the pion form factor [isf. Together, the two appear to support the pQCD predic- 
tion of as/IQ^I variation of the form factor at large \Q'^\- Bebek et ai.0| have reported 
|Q2|F^(9.77 GeV^) = 0.69±0.19 GeV^ for the spacelike form factor. Within errors this is 
consistent with being nearly factor two smaller than the timelike form factors for > 9 
GeV^, as found for protons. 

Our measurement of the kaon form factor stands alone at present. The asymptotic pQCD 
prediction is [5] F^^kW\) = Svra J2^_^+/|Q2|. Por a, = 0.3, = 130.7 ±0.4 MeV, and 
fK+ = 159.8 ± 1.5 MeV O, this leads to F^(13.48 GeV^) = 0.010 and Fi^(13.48 GeV^) 
= 0.014, which are factors of ~8 and ~4 smaller than our results, respectively. The 
and Q2 independent pQCD prediction F^{\Q'^\) / Fk{\Q'^\) = fl/fl = 0.67±0.01 is also in 
disagreement with our result F^(13.48 GeV2)/Fi^( 13.48 GeV^) = 1.19±0.17. Bakulev et 
ai.[l5| have estimated soft contributions to the pion form factor in the framework of QCD 
sum rules. Their formulation leads to F^(13.48 GeV^) = 0.010 and F^{\Q^\)/ Fk{\Q'^\) = 
0.51, i.e., the discrepancy between our experimental results and the theoretical predictions 
remains. We note that this behavior is in contrast to the good agreement between the 
measured tt^'J transition form factor and the pQCD prediction for the same jl6l |. 

Our result for (7^^(13.48 GeV^) is in excellent agreement with the results of the Fermilab 
E760/E835 experiments in which the reverse reaction pp — > e+e~ was measured j^. Our 
results provide the first independent confirmation of the Fermilab results. 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. This work was supported by the National Science Foun- 
dation and the U.S. Department of Energy. 
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